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Abstract 

Verhulst logistic curve either grows OR decays, depending on the growth rate parameter value. 
A similar situation is found in the Gompertz law about human mortality. However, growth 
can neither be infinite nor reach a finite steady state at an infinite asymptotic time. Moreover 
before some decay, some growth must have occurred. For example, some levelling-off could 
occur at finite time, followed either by some growth again or then by some decay. Numerous 



examples show that the Verhulst and Gompertz modelisations are too reductive (or restrictive) 
descriptions of their original purpose. It is aimed, in the present note, to encompass into ONE 
simple differential equation the growth AND decay features of, e.g., population sizes, or numbers, 
but also of many other measured characteristics found in social and physical science systems. 

cc 

Previous generalisations of Verhulst or Gompertz functions are recalled. It is shown that drastic 

i growth or decay jumps or turnovers can be readily described through drastic changes in values 

i— I of the growth or decay rate. However smoother descriptions can be found if the growth or decay 

> 

ON rate is modified in order to take into account some time or size dependence. Similar arguments 

can be carried through, but not so easily, for the so called carrying capacity, indeed leading to 

ON more elaborate algebraic work. 

O 
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I. INTRODUCTION 



Verhulst 1845 seminal paper pQ bore on previous (1798) considerations by Malthus [2] on 
population size evolution. Verhulst model attempts to describe how the. growth (or mutatis 
mutandis its decay) of a single species towards its equilibrium population is limited, taking 
into account the reproduction rate r and the nowadays so called limiting capacity M, i.e., 
the maximum value the population can reach in its environment., due to limited resources. 
He obtained what he called the logistic curve (or map). In a subsequent paper, Verhulst [3] 
attempted a different mathematical form of the blocking growth term, thereby obtaining an 
exponentially slow approach of the carrying capacity, but not a sigmoid. 

There have been many considerations starting from Verhulst approach of population size 
evolution, generalising his original work. However, it is of interest to recall that among the 
generalisations, one should distinguish between those which stick to one single differential 
equation, and those which involve a set of equations, because they consider a competition 
effect with another population. Among the former ones, I am fond of the 1939 work by 
Avrami j3] on crystal growth, - the so called species changing its thermodynamic phase 
or state, but there are many others. Among the latter, one has only to go back to the 
work of Lotka, published [5], in 1925, and of Volterra [6], in 1931, pertaining to so called 
prey-predator models. Not attempting to cite here many subsequent publications., let it be 
pointed out that Blanco [7] has suggested a trivial relationship between the logistic equation 
and the Lotka- Volterra (IV) models, indicating that the former is in fact a sort of limit of 
the LV equations, pending a redefinition of the parameters. 

However, growth can neither be infinite nor reach a finite steady state at an infinite 
asymptotic time, as in the logistic map. Moreover before some decay, some growth must 
have occurred. In particular, some levelling-off could occur at finite time, though followed 
either by some growth again or then by some decay, in both cases, with a smooth or an 
oscillatory variation. 



It is of interest in contrast, to recall ideas byGompertz [8] who in a series of papers 
near 1825 discussed 'the Nature of the Function Expressive of the Law of Human Mortality, 
and (...) a New Mode of Determining the Value of Life Contingencies". The Gompertz 
model for describing the decreasing number of members of a population can be in fact easily 
reformulated in terms of a (simple first order differential) growth equation of Malthusian 
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time 



FIG. 1: Comparing a Verhulst and a Gompertz growth or decay law (lhs), and their first derivative, 
in the case of parameter r > for growth (rhs); observe the asymmetry in Gompertz laws 

form, with an exponentially decaying birth rate. This leads to another form of an analytical 
sigmoid curve, though asymmetric, in contrast to the Verhulst logistic map. 

In brief, depending on the sign of the parameters in Verhulst or Gompertz approach, 
one can describe some growth or some decay see Fig. [T] (lhs). However numerous examples 
show that both Verhulst and Gompertz modelisations are thus too reductive (or restrictive) 
descriptions of their original purpose. Whence, it is of interest to aim at finding some simple 
description combining some less trivial behaviour, i.e. obviously, practically inducing some 
extremum at some time in the function which should be the solution of a, by request, simple 
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first order differential equation. This is the aim of this note, i.e. to encompass into ONE 
simple first order differential equation growth AND decay features of, e.g., population sizes, 
or numbers, but also of many other measured characteristics found in social and physical 
science systems. 

In Sect. 2, after recalling Malthus and Verhulst equations, a few applications are men- 
tioned in order to substantiate the broad framework in which Verhulst equation can be 
applied. Growth, indeed, is not only found in population studies, through counting the 
number of individuals in some system, but can be measured through other means, like the 
number of stems, the mass, the length, ... of "something". It is often found, of course, that 
the Verhulst sigmoid is only a rough representation of the data. Whence one has sometimes 
invented artificially skewed logistic curves j9]; see Sect. 2, for a general introduction to 
growth and such subsequent considerations. 

Subsequent considerations on decay, with examples about languages and religions, are 
to be found in Sect. 3, before outlining Gompertz law. Discussion about bio-demography 
applications can be found in an interesting 1997 review, up by Olshansky and Carnes [10J. 
Nevertheless, ad hoc, previous generalisations of Verhulst or Gompertz functions are briefly 
recalled to conclude Sect. 3. 

In Sect. 4, a few examples show that the Verhulst and Gompertz modelizations are indeed 
too reductive (or restrictive) descriptions of reality. However, drastic growth or decay jumps 
or evolution turnovers can be readily described through drastic changes in values of the 
growth or decay rate, as shown inSect. 5. Sometimes the evolution show drastic jumps 
or drops indeed, - due to exogenous or endogenous causes. However smoother descriptions 
must be found, i.e. combining some less trivial behaviour. Obviously, one should aim at 
practically inducing some extremum at some time in a smooth function which should be the 
solution of a, by request, simple first order differential equation. 

It is aimed, in Sect. 6, to encompass into such ONE simple differential equation the 
growth AND decay features of population sizes, or numbers, but also of many other measured 
characteristics found in social and physical science systems. In fact one can observe that one 
could consider that he growth or decay rate has either some time or some size dependence. 
The algebra is rather simple in order to recover a Verhulst or Gompertz evolution equation.. 
In Sect. 7, the time and/or size dependence of the so called carrying capacity is considered 
in the same spirit as the growth rate. However this is less trivial, indeed leading to more 
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elaborate algebraic work. Sect. 8 serves as a conclusion. 



II. INTRODUCTION TO GROWTH 

Following the above historical points, let us do some mathematics and data analysis 
within Verhulst equation framework. 

A. Growth equations 

Recall Malthus model or rather his equation for describing a (normalized) population 
size x evolution as a function of time t, due to a birth rate r, i.e. 

- = rx, (1) 

which has the (catastrophic) exponential solution x(t) = e rt . 

The original Verhulst population growth model considers that resources are limited, such 
that the evolution of the number of members N(t) in the population is described in terms 
of a carrying capacity M and a growth rate r, according to 

N 



dN 

— = r N 
dt 



(2) 



With the change of variable x = N/M, if M is a constant, a dimensionless form equation 
can be thereby written, i.e. 

— = r x [1 — x] . (3) 
It has for solution, the so called logistic map, a sigmoid curve, 

rt 

* = T^— t > ( 4 ) 

if r is a constant as well. Notice that Verhulst had examined complicated forms of the 
growth limiting term, like a square, but also a square root [T], both cases without any theo- 
retical justification, but aiming at better fitting available data Indeed, there is no biological 
justification for this assumption; Gilpin, Case, & Ayala [TT] proposed a growth equation of 
the form dx/dt = r x[l — (x/M) e , - a larger than 1 value of 9 resulting in a steepening of 
the growth curve and a higher value for the inflection point, as easily calculated. 
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B. Applications of the logistic map 

Let us point out to a few cases of the "Verhulst population growth model", as discussed 
during the oral presentation, in order to show some universality value of the logistic map 
application. 

• Population size: Beside those considered by Verhulst, an interesting application on 
population size is that of Montroll and Badger [12], in 1974, who studied the USA 
population between 1810 and 1970. An interesting parameter value is found, i.e. 
r = 0.02984, which implies interesting considerations on american sexuality during 
those times. 

• Mass growth: However the logistic map can also be applied to, e.g., mass growth, 
like in the report of Lee and Werning [T3j, in 2008, who compared Alligators, Ele- 
phants, Shrews and Ostriches, on one hand, but in the other hands Tenontosauruses, 
Allosauruses, and Tyrannosauruses, using skeleto chronology; see reference [13] for the 
parameter r values, in kg/yr, and subsequent discussions. For completeness, let us 
mention that the growth curves were compared with an alternative model to that of 
Verhulst, i.e. the von Bertalanffy [T4] curve of mass accumulation. 

• Length growth: As far as 1928, Pearl group studied the Cucumis Melo stem length, 
either under field conditions [15] or in the absence of exogenous food and light [T6] . 
The population was considered to be the number of nodes. Notice that r is of the 
order of 0.7 in their observations. 

C. Skewed logistic 

In so doing, i.e., studying Canteloup growth, Pearl et al. empirically found that one 
should generalised the logistic map in order to reproduce better the data. Rather than the 
mere logistic, they propose a form like 

k 

V® = 1 _|_ e a +a x t+a 2 t*+... ' ( 5 ) 

where y, in [16], is the number of seedlings of the canteloups. 
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It is quasi obvious that a more general form is like 



y(t) 



k 



(6) 



1 + e/W 



to be discussed below. 

III. INTRODUCTION TO DECAY 

Recent applications along Verhulst (limited growth) approach have been concerned with 
the decay of a population size, rather than its growth. 

A. Religions 

The time evolution of several "main" religions was considered to be described, at a so 
called microscopic level, along the lines of the Avrami-Kolmogorov equation, describing a 
liquid-solid state formation in a continuous time framework [17], i.e. a peculiar form of 
Verhulst law [18J. The solution is usually written as 



where c(t) is the volume fraction being transformed from one phase to another; K and n are 
adjustable parameters. However, it is "obvious" that before decaying the number of adepts 
of any religion had to grow from the founder, N(0) = 1. Nowadays, several "religions" have 
still a growing number of adepts. Others have reached a maximum, thus decay, or go again 
into a growing pattern. Population growth-death equations have thus been conjectured for 
modelling such religion size evolution dynamics, in a continuous time framework [T9] . [20J, 
[2"T] . [22], [23] and [24J attempting to describe the existence of maxima in the number of 
adepts. An interpretation of the features, after data fitting, can resort to an attachment 
parameter, as resulting from joining at first one's "mother religion", within exogenous field 
conditions, as one "joins" at first one's mother tongue. 

B. Language death 

For completeness, notice that the disappearance of world's languages has also been re- 
cently considered through a Verhulst time evolution equation in terms of the number density 



c(t) = 1 - e 



-Kf 



(7) 
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Pi of practitioners of some language i, in [25] and [26], i.e., 

dpi f, Pi 



at r «H 1 -5U (8) 

where C{ is the carrying capacity of the environment for the population speaking the i 
language and is a "negative growth" rate. The values are at this level of no interest. 
However notice that the speakers do not fully represent a single so called population, since 
their member evolution is toward a bilingual attitude. In fact, there is a competition between 
two (or more) languages/populations, - which according to the point of view mentioned in 
the Introduction rather pertains to a prey-predator problem. 



IV. GOMPERTZ DEATH LAW 

In fact, previous to Verhulst, Gompertz [5] discussed human population mortality dynam- 
ics from a practicing actuary point of view. Gompertz used equal long time intervals, i.e. 
10 year long intervals, to observe that the differences between the common logarithm of the 
number of persons living in such successive equal age intervals were almost identical during 
a significant portion of their life span. Thus the numbers of living individuals in each yearly 
increase of age are in geometrical progression. Whence mortality progresses geometrically 
as age increases arithmetically. Therefore Gompertz proposed a 3 parameter expression for 
the number N(t) of survivors of a population at time t: 

N(t)/N inf = exp[-b exp(-rt)}, (9) 

where N in f, b and r are positive constants. This corresponds to an exponentially decaying 
birth rate r in Malthus equation, i.e. r = r e~ Kt , pending r and k being positive constants. 
For completeness, notice that Makeham introduced an age independent mortality rate to be 
considered in order to better represent available data [27J, i.e. adding a A + Bt term in the 
right hand side of Eq.(|9j thus attempting to take into account extrinsic effects. 

The easily derived differential equation describing the exponential rise in death rates 
between sexual maturity and old age, is commonly referred to as the Gompertz equation, 
i.e., 

rfnr h 

(10) 



dx , 
— = r x Loq 
dt 
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where k has mutatis mutandis the same meaning as M, the carrying capacity, in the Verhulst 
approach. Notice that Gompertz and Makeham wondered whether k can be interpreted in 
a socio-biological framework, i.e., as a biological limit to life, based on some vital force. 



A. Skewed Gompertz law 

In the original Gompertz model, r and k are supposed to be constant; b, in Eq.((9|), and k, 



in Eq.( 10 ), are related to each other through r and the initial (or final) conditions on N(t). In 
1924, Trachtenberg [28J already examined extensions of Gompertz law, sometimes departing 
from Gompertz geometric progression hypothesis, - extensions which do not necessarily lead 
to a sigmoid curve, but can be called skewed Gompertz laws. One of them reads 

N(t) = N inf e l"+M+<* 2 +dt 3 +et*] (U) 

It can clearly be the solution of a second order differential equation for log(N(t)), as wished 
by Gompertz, and presents two inflexion points at precise t a and tp. 

In fact, similarly to Eq.([6]), the most general extension of a Gompertz form can be hereby 
proposed to be 

N(t) = N inf exp[-b exp(-g(t))] . (12) 



B. Biomedical considerations 



For the purpose of the reader general information, an interesting 1997 review by 01- 
shansky and Carnes [TO], on "ever since Gompertz", should be quoted. However, since the 
discussion is mainly about demography applications and biological causes of mortality, in 
various species, - like the vital force concept and in particular the influence of radiations, the 
paper content is quite aside our present considerations, and not further elaborated upon, 
though of great interest. 

Nevertheless, it is of common knowledge that exogenous causes can manifest themselves 
and influence some "decay" in the number of living individuals in a population. Take the 
trout population data of Lake Ontario measured by Beeton [29], quoted by Meadows et 
al. [30J, for example. Pollution is supposed to be the influential external field. On the 
other hand, a case of endogenous-like cause was noticed by Kormandy [31], in 1969, who 
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commented on data, which I did not find, and who stated that deer or goats, even though 
natural enemies are absent, "often overgraze their range and cause erosion or destruction of 
the vegetation", and consequently die. Finally, on the Bauduen, Var city wikipedia website 
{http : // fr.wikipedia.org/wiki/Bauduen) it is reported that the city population, since 
1793, increased till 1820, thereafter decreased till 1970 and is going back up nowadays. In 
this case, sharp turnovers are found at the time of war ends. 

V. GROWTH AND DECAY AND GROWTH AND DECAY AND ... 

In the case of Lake Ontario trouts, here above mentioned [29] , there is obviously a max- 
imum near 1925, followed by a rather sharp drop of three orders of magnitude in 1955 or 
so. Similar though smoother, features are found for Bauduen population. In brief, such few 
examples should remind us that, in many measurements, when counting the evolution of a 
population (of whatever), there is a succession of maxima and minima, though the growth 
and decay regimes can be of different sizes and time spans. Remember business cycles as 
well [32], [33], [31], [35], or [36]. To debate on the evolution laws per se is out of question 
here, but to find an equation containing both types of extrema is the goal of this report, 
basing our considerations on either Verhulst or Gompertz equation. 

A. Beyond Verhulst and Gompertz equations 

In view of the above, it is obvious that many multiparameter extensions have been pro- 
posed in order to "improve" Verhulst or Gompertz laws [37], see App. A. However, parameter 
justification is often missing, within our present requirements. Thus a simple mathematical 
framework to find several inflexion points on a time dependent "growth" curve is still of 
interest. One may imagine to introduce extra terms, representing endogenous or/and exoge- 
nous causes for growth or decay in Verhulst or Gompertz formalisms, Let us mention Bass 
model [3B] of innovative and imitative behaviour of consumers, or Zwanzig introduction of 
a (damped exponential) heredity kernel [31] , both generalising Verhulst approach. There is 
to my knowledge no equivalent extension within Gompertz framework. However both, Bass 
and Zwanzig considerations, have more physical bases than other purely mathematically 
oriented extensions. 
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For completeness, let us recall Bass model [38J of new product acquisition by consumers; 
it leads to a single "diffusion" equation with constant coefficients, - perfectly suitable to our 
aims: 

d * = [N-X(t)}[p+(^)X(t)]. (13) 

where X(t) is the number of product users at time t; N is a number of potential users, p is 
the coefficient of innovation, q is the coefficient of imitation. It is a Verhulst-like equation, 
though the r.h.s has an extra zero order term. The solution reads 

1 _ e -(p+g)t 

X(t) = N — - — -. — - . (14) 

By differentiation, one can easily find the time at which the function has a maximum and 
its value. 

On the other hand, Zwanzig [35] writing of Verhulst equation amounts to a change of 
variable e v = N/M, - somewhat as in Montroll [ID] considering that the evolution concerns 
log(N) rather than N, and the introduction of a heredity kernel K{t) = K(0) e t//T , such 
that one obtains an integro-differential equation (not written here). K(0)T is in fact like 
the growth rate r in Eq.Q. For finite T, one obtains an oscillatory but damped evolution, 
with an overshoot of the carrying capacity value in fact. 

Both previous formalisms present a hump in the "population" evolution. Nevertheless 
already within e.g. Verhulst framework, it can be observed that describing either (i) sharp 
turn-overs or (ii) drop or jump transitions is possible. It is for example, simple to introduce 
a drastic change in the r value. Similar cases can be illustrated with Gompertz law, as 
shown here below. 



B. Sharp turn-overs 



The four possible cases of sharp turn-overs forcing a continuous evolution of the popu- 
lation number, in a Verhulst approach, due to a drastic shift in r values, i.e. from some r 
to some — r value, or the opposite, at some time t are shown in Fig. [2j Similar curves are 
shown for the Gompertz evolution law on Fig. |3j Notice the pronounced asymmetry in the 
latter case. 
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time time 



FIG. 2: The four possible cases of sharp turn overs with continuous evolution in a Verhulst 
approach due to a drastic shift in r values at different times t 

C. Drop or jump transitions 

On the other hand, the four possible cases of population drops or jumps in a Verhulst 
approach due to a drastic shift in r values, at some fixed time t are shown in Fig. |4j as 
if a strong external event had influenced a shift from some r to some — r, or the contrary. 
Similar curves are shown for the Gompertz evolution law on Fig. [5j 
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FIG. 3: The four possible cases of sharp turn overs with continuous evolution in a Gompertz 
approach due to a drastic shift in r values at various times t 

D. Time and size roles 



However such drastic change of signs of the growth rate are rather unlikely or very rough 
approximation of reality. A smoother evolution of r must be searched for. In this respect, 
let us recall that the dimension-prone Verhulst differential equation, Eq.Q, can be easily 
integrated, 
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FIG. 4: The four possible cases of population drops or jumps in a Verhulst approach due to a 
drastic shift in r values at some fixed time t (=+1 , or -1) 

• if the growth rate r is constant, 

• and if the carrying capacity M is a constant. 

Let us first observe that there are two equivalent ways to write both Verhulst and Gom- 
pertz equations, either emphasising the time role or the size role. 
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FIG. 5: The four possible cases of population drops or jumps in a Gompertz approach due to a 
drastic shift in r values at some fixed time t (=+1 , or -1) 

I. Time role 



It is easy to observe that, after an appropriate simple change of variables, the correspond- 
ing differential equation to Eq.Q, e.g. for the skewed logistic, Eq.(j6]), could read 



dy b 



y [k-y] v(t) 



(15) 
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- thereby defining v(t). One can easily integrate both sides of the equation, by the technique 
of separation of variables, i.e., 

dy b 



y [k-y] k 

and easily obtain the skewed logistic therefore. 



v(t) dt , (16) 



2. Size role 



Another interesting way, for further discussion, is to rewrite the evolution equation, letting 
the time t in the right hand side be an implicit variable, thus referring to the evolution 
through the size itself, i. e. let 

t - v ^ < 17 > 

where V(y) = ry(l — y) in Verhulst model. Imposing the stability condition dV Jf^ = V(y) = 

0, one easily finds out that the r.h.s. is like 

V(y) = ay + by 2 + cy 3 + (18) 

1. e. there is no zero order term. In other words, keeping the cubic term only, leads to 

dy 

— = ay(y - y b )(y - y c ) , (19) 

therefore simply generalising the Verhulst equation leading to the logistic map, thereby 
introducing two inflexion points at y b and y c . 



E. Second order derivative 



Without going into elementary considerations for mathematicians and physicists about 
the conditions on a differential equation to produce, as its solution, a curve with an ex- 
tremum, it is immediately obvious, in the present framework that to differentiate once more 
Eq.Q and Eq.(j9| will produce a maximum, because of the introduction of a second inflexion 
point; - see Fig. [l] (rhs). However, the goal is to find a first order differential equation 
(ODE) not a second ODE. As hinted here above, a time or size dependence of r and/or k 
(~ M) seems the simplest way. 
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VI. GROWTH RATE ROLE 



Growth rate role in Verhulst model 



One may rewrite Verhulst equation, Eq.([3]), explicitly expressing the time dependence on 
the right hand side, i.e. 

-rt 



dx 
~dt 



re 



1 + e~ rt ) 2 2(1 + coshirt)) ' ^ 
In order to have a curve with an extremum the most simple way seems to consider its 
derivative; the extremum being located at the inflexion point of the logistic map. One finds 



d cc r _ _ dec 
— = r 1 - 2x] — 
dt 2 1 J dt 



r 2 x (1 — x) (1 — 2x), 



or emphasizing the time dependence on the right hand side, i.e. explicitly writing 

d 2 x r 2 e rt (1 - e rt ) 



dt 2 

One can also rewrite the above as 

d 2 x _ r 1 - ~ rt 
dt? " 



rt\3 



1 + C 



x (1 — x) = Rv x (1 — x) 



(21) 



(22) 



(23) 



1 + e rt _ 

in order to reproduce the r.h.s. of Verhulst equation, but with a new r, called Ry having 
a specific time dependence, i.e. as 

,1 -e rt ~ 



R\ 



—r tanh{rt) 



(24) 



1 + e rt 

- which is not a birth-death growth rate of any known population, up to now. In so doing 
one obtains the "generalized" time- Verhulst equation, admitting that the growth rate should 
be time dependent. 

Considering a writing in terms of a size dependent r, one can write 



d 2 x 

du 



d , dx ^ 
dt^~dt' 



Pv 



,dx, 



i 1 - & 



which has the Verhulst form ^ = py u [1 — u], with, comparing Eq.(21) and Eq.(23), 



Pv 



Rv/r 



r (1 - 2x) 



(25) 



(26) 



1 — r x (1 — x) 1 — rx(l — x) 
where pv is expressed in terms of the population size. If one wishes py to be expressed in 
terms of the time t, one should replace u (= x) and Ry by their expression found in Eq.Q 



and Eq.(24) respectively; e. g., Ry — r (1 — 2x{t)). 
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B. Growth rate role in Gompertz model 



Similar manipulations and algebra can be performed for the basic Gompertz model, i.e., 



log 


~k 




x log 


k 




- 1 






X 






X 



d 2 x 
df? = 

thereby defining R G , in terms of its size dependence, 

Rg 

Therefore, in order to write 



R G x log 





~k 




log 




- 1 




X 





d 2 x 



d dx s 
d~V~dl' 



Pa (f ) log 



k 



(dx j dt) 

which has the Gompertz form, one has to introduce p G with 

Rc/r log [k/x] — 1 



Pg 



(27) 



(28) 



(29) 



(30) 



log [k/x] — log [r log(k/x)] log [k/x] — log [r log(k/x)] 
where pc is expressed in terms of the population size. If one wishes pa to be expressed in 



terms of the time t, one should replace Rg and (dx/dt) by their expression given in Eq.(28) 
and Eq.([9]). It seems that such a writing is out of necessity here. 



VII. CARRYING CAPACITY ROLE 



Beside the possibility of expressing the growth rate in terms of some time or size de- 
pendence, let us recall that k and M in Gompertz and Verhulst model respectively are 
considered to be constant. To remove such a restriction has already been considered in 
works referenced to here above. However such a consideration is of highly relevant modern 
interest in our growing world, and thereby requests some attention. 

Notice that some population growing in a limited environment can approach the ultimate 
carrying capacity of that environment in several possible ways: (i) it can adjust smoothly to 
an equilibrium, (ii) it can overshoot the limit and then dies back again either in a smooth 
very damped way or (iii) in an oscillatory, more or less damped, way. Let it be emphasised 
that the population can overshoot the carrying capacity and in the process decreases the 
ultimate carrying capacity by consuming some necessary nonrenewable resource: e.g., deer 
or goats, when natural enemies are absent, often overgraze their range and cause erosion or 
destruction of the vegetation [31] 
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FIG. 6: Possible cases of population evolution in a (lhs) Verhulst or a (rhs) Gompertz approach 
when the carrying capacity has a time dependent component, as compared to an ordinary growth 
or decay evolution for constant r or k (~ M); r = 1 for all curves; used analytic forms are given in 
the main text. A bump or a dip can be seen to occur here at negative finite time 

A. World model 

Related considerations on the carrying capacity are used in Meadows et al. [30] on the 
"standard" world model. At first it assumes no major change in the physical, economic, 
or social relationships that have historically governed the development of the world system. 
Food, industrial output, and population grow exponentially, as according to Malthus [2J, 
until the rapidly diminishing resource base forces a slowdown in industrial growth. Popula- 
tion growth is finally halted by a rise in the death rate due to decreased food and medical 
services. However, because of natural delays in the system, the population (and pollution) 
continues to increase for some time after the peak of industrialisation. 

Developments of the world model have considered the doubling or resources [30], - to- 
gether with plenty of other interfering causes, for describing the population evolution. 
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B. Carrying capacity role in Verhulst and in Gompertz model 



The time and/or size dependence of the so called carrying capacity could be considered 
in the same spirit as the growth rate. However this is less trivial and apparently can only be 
tackled through numerical work. Indeed, one can separate the variables when one has either 



a form k(x) or k(t) in Gompertz model, or M(N) or M(t) in Verhulst model, in Eq.(lO) 



and Eq.(15) respectively, i.e. 

dx , s , \k(x(t)) 
— = r x(t) log 
dt w y 

dN r 



(31) 



N(t) [M(N(t)) - N(t)] , (32) 



dt M(N(t)) 

However to go on, some analytic form for k and M is needed. This is leading to more 
elaborate, less universal, algebraic work, and is outside our present report. 

Nevertheless in order to substantiate the role of the carrying capacity time dependence 
effect in the Verhulst and in Gompertz model, a posteriori analytic cases have been consid- 
ered. A few results are shown in Fig. [6] when a time dependent behaviour is a posteriori 
included as a mere factor into the original Verhulst of Gompertz laws, r being kept constant. 
For completeness, the used equations in Fig. [6] read 

/ x 1 r e* r e* - 1 
Mt) = - 2 [j^] + 1] (33) 

= lie-'") [±4^| + 1] ■ (34) 

2 e l + 1 

Bumps or dips are markedly seen, within overall growth or decay features, even in this 
simple manipulation. A size dependent effect of the carrying capacity either in the Verhulst 
or Gompertz approach has not been studied. 



VIII. CONCLUSION 



A goal has been to describe some simple way toward one single simple first order dif- 
ferential dynamics equation, recalling that 3-parameter families of curves often describe 
growth data adequately, like in Verhulst or Gompertz approach. It has been stressed that 
extensions of the corresponding equations in order to introduce extrema, or more than one 
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inflexion point, are easily found in the second order derivative. However a change of vari- 
ables is necessary, implying that the growth rate or/and the carrying capacity becomes time 
(or size) dependent. This corresponds in practice at considering exogenous or/and endoge- 
nous causes, in other words, external or internal field interactions with the population, i.e. 
intra-community interactions in the latter case. 

Bass model [38] and Zwanzig model [39] already contain such extrema in an ad hoc way 
within Verhulst population growth framework. Here it has been aimed at providing a single 
diffusion Riccati equation with time/size dependent coefficients of realistic nature, enforcing 
that the corresponding first ODE is of the Gompertz or Verhulst type. An analytic form 
for the growth rate coefficients has been derived. The carrying capacity form seems not so 
immediately tractable and demands more numerical work. 

Other ways of interest to describe a population evolution are those taking into account 
some time delay, as in the world model, and in work on anticipation by Dubois [JT] or 
Rosen [12] . Indeed, in biological and sociological terms, the logistic and Gompertz equations 
present a logical absurdity, since it requires that cause and effect occur quasi simultaneously. 

In fine, let it be pointed out that considerations of discrete time effects do not seem 
to have been studied on Gompertz model: that might be of interest in the framework of 
CASYS conference reports. 
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Appendix A. Multiparameter extensions 

For completeness, it seems appropriate to recall generalisations of Verhulst and Gompertz 
laws, beside the so called skewed laws, mentioned in the main text. E.g., one can write indeed 



where p, a, and 9 are to be positive numbers, - the logistic map corresponding to 6=1. 




= e 



-pe 



(35) 
(36) 
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Ahuja and Nash [37] have further generalised these distributions by introducing an addi- 
tional parameter cf) > 0. The density functions of these distributions are respectively 



However, it seems that there is no physical justification of these parameters to be usefully 
discussed within the present context. 
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